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Abstract 

We present an overview of the options for diffraction implemented in the general- 
purpose event generator Pythia 8. We review the existing model for low- and 
high-mass soft diffraction and present a new model for hard diffraction in pp and 
pp collisions. Both models uses the Pomeron approach pioneered by Ingelman and 
Schlein, factorising the single diffractive cross section into a Pomeron flux and a 
Pomeron PDF. The model for hard diffraction is implemented as a part of the 
multiparton interactions framework, thereby introducing a dynamical rapidity gap 
survival probability that explicitly breaks factorisation. 


Presented at the 16th conference on Elastic and Diffractive scattering (EDS Blois 2015) 



1 Introduction 


While most phenomena in high-energy hadronic collisions have been explained by QCD, the effects 
of the softer hadronic collisions remains a mystery. We observe these collisions in experiments, 
and can motivate why they should be present, but the explanation of how they occur is still largely 
based on phenomenological models. These models should be able to describe all aspects of such 
collisions, like differential cross sections, one-particle distributions and global event characteristics. 
The models should also describe the exclusive topologies of these softer collisions, specifically the 
occurrence of rapidity gaps. 

Many models, including the models used in Pythia 8 [I], are based on Regge theory. In 
this theory, poles in the plane of complex spin a can be seen as hadronic resonances. These 
appear to lie on linear trajectories, a{t) = a(0) + a't. Most important for high-energy collisions 
is the Pomeron (P) trajectory, with its a(0) > 1 explaining the rise of the total cross section. 
This state is a colour-singlet carrying the quantum numbers of the vacuum. From a modern 
viewpoint it (predominantly) consists of gluons and could thus be called a glueball, or a gluonic 
ladder if present in the final state (a cut Pomeron). A topological expansion can be defined, 
with increasingly complex processes. The simplest possible exchange is a single-Pomeron one, 
which gives rise to elastic scattering. Multi-Pomeron exchange is also possible, e.g. involving the 
triple-Pomeron vertex. This way various diffractive topologies can be constructed. In this paper 
we focus on the single diffractive (SD) topologies, since these have the largest diffractive cross 
section and form the starting point on the road towards more complex configurations. 

Ingelman and Schlein [2] proposed a model in which the exchanged Pomeron can be viewed as 
a hadronic state. This opened up the possibility for using Pomeron parton distribution functions 
(PDFs) to be combined with a probability for taking out a Pomeron from the initial hadronic state, 
the Pomeron flux. The diffractive system can be viewed as a hadron-hadron collision at reduced 
energy, and existing hadron-hadron event generators can be used for modelling the diffractive 
events. The simplest model does not allow for multiparton interactions (MPIs), however, or 
equivantly for the final-state effects of multiple cut Pomerons. These MPIs create additional 
colour strings in the event, each string giving rise to hadronic production. Hence we risk filling 
up the rapidity gap created by the exchange of the ‘first’ Pomeron. As a rapidity gap is needed 
to trigger on diffractive events, we risk losing a large fraction of the could-have-been diffractive 
events by these MPIs. This introduced the concept of rapidity gap survival probability (RGSP), 
which is unique to hadron-hadron collisions, given credibility by the lower observed rate of hard 
diffractive processes at the Tevatron than expected from HERA flux/PDF determinations [3]. 

2 Soft diffraction in Pythia 8 

The soft diffraction machinery available in Pythia 8 was originally developed for Pythia 6 
[1] , but rewritten and expanded for the new version, and now includes both single-, double- and 
central-diffractive systems (SD, DD, CD) as well as elastic collisions and non-diffractive topologies 

[5] . The total hadronic cross section is calculated using the Donnachie-Landshoff parametrisation 

[6] , with a Pomeron and Reggeon term. The elastic and diffractive cross sections are based on 
the Schuler-Sjostrand model [7j and the non-diffractive cross section is inferred from these two 
models. 

The Schuler-Sjostrand model is also based on Regge theory and gives an approximate dM/M^ 
mass dependence as well as an exponential t dependence. Fudge factors have been introduced to 
the model, to dampen the cross sections close to the kinematical limits, as well and to dampen 
the DD cross section where two diffractive systems overlap. Other Pomeron-flux models have 
also been implemented in Pythia 8 (see the manual [8]). The subsequent hadronisation of a 
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diffractive system is a separate chapter, and the same for all Pomeron-flux models. Particle 
production depends strongly on the mass of the diffractive system, however, and hence it has 
been split into two regions. 

2.1 Low—mass soft diffraction 

In the low-mass regime, M < 10 GeV, energies are not sufficiently high to apply a perturbative 
framework to the Pomeron-proton subcollision. Instead we visualise the event as an interaction 
where the Pomeron has “kicked out” a parton from the diffractively excited hadron. If a valence 
quark is kicked out then a single colour string is stretched between it and the diquark remnant. A 
kicked-out gluon gives a hairpin string topology, stretching from one quark in the proton remnant 
to the gluon and then back to the remaining diquark of the remnant. The probability for the 
Pomeron to interact with either a quark or a gluon is mass-dependent, P{q)/P{g) = N/M^ with 
p being a tunable parameter, making the gluons dominate at higher mass. There are no additional 
MPIs in the low-mass regime. The strings are hadronised using the Lund string fragmentation 
model [9] and gives rise to low-pT activity in the diffractive system. 

2.2 High—mass soft diffraction 

In the high-mass regime, M > 10 GeV, a perturbative description is attempted. So as not to 
give any discontinuies, and possibly also representing a real physics evolution, the fraction of 
perturbative events gradually increases with M and dominates for M > 20 GeV. 

In the new component the Pomeron is viewed as a particle with partonic content a la Ingelman 
and Schlein. Thus, once M and t have been selected, the system is set up as a Pp collision and 
a semi-hard perturbative 2 —>• 2 partonic interaction is selected by the MPI machinery. Inside 
the Pomeron-hadron system the full interleaved evolution of initial- and final-state showers (ISR 
and FSR) and MPI is applied using the Pomeron PDFs. The MPI activity in the subsystem has 
been tuned to give approximately the same amount of activity as in non-diffractive events of the 
same mass, by introducing an effective total Pomeron-proton cross section. This (tunable) total 
cross section is set to a constant value of 10 mb, slightly higher than other numbers found in the 
literature. The colour strings obtained in the evolved diffractive system are hadronised using the 
Lund string fragmentation model. Jets can be produced in the 2—7-2 partonic processes. 

Although the models for soft diffraction available in Pythia 8 are largely successful, some 
minor issues show up. Not all aspects of the data are described using the default model and 
settings, both on the level of differential cross sections and on that of particle spectra. A retune 
of parameters used in the default model could fix some issues, in particular if allowing for more 
flexible shapes e.g. for the Pomeron flux. We intend to improve the default models in the near 
future. 

3 The new model for hard diffraction in Pythia 8 

The model described above does allow for QGD 2 —)• 2 processes at all pT scales, but is primarily 
intended for lower pT values. It is not intended for the study of truly hard processes, either 
in QGD or beyond. Instead a model for hard diffraction has been developed [10] based on the 
assumption that the proton PDF can be separated into a non-diffractive and a diffractive part. 
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with the diffractive part described using the factorisation approach, 


Q^) = 


1 rt „ 



dt dx' dxp /p/p(xp,t) /i/p(x',Q^) 5{x - x'xp) 


= / — /p/p(a;p) /i/p( — 
Jx Xf> xp 


( 1 ) 


The probability of diffraction on one side is then given as the ratio of diffractive to inclusive 
PDFs, 

pD(x„ g2) = Q2 )//,/p(x„ Q2). (2) 

At high energies most interactions occur at low x where P^(x,Q^) ~ 0.1. Hence we expect 
approximately 10-15% of the events to be diffractive based only on Eq. 

In addition the model implements a dynamical gap survival. This means we do not allow any 
further MPIs to occur between the two incoming hadrons, so as to ensure the gap survives. In 
practise the tentative classification as diffractive, based on Eq. initially has no consequences: 
all events are handled as non-diffractive hadron-hadron collisions. Only if no additional MPIs 
occur does a diffractive classification survive and only then is the Pp subsystem set up. A full 
evolution of ISR, ESR and MPIs is performed in this Pp system, along with hadronisation of the 
colour strings in the event. At this stage all non-diffractive events can be discarded for a pure 
diffractive sample, or can be kept and hadronised as usual for an inclusive sample. 

The restriction on the number of MPIs in the hadron-hadron system introduces an additional 
suppression factor of ~ 0.2. With this method we can explain the observed “factorisation breaking” 
at the Tevatron, without introducing any new parameters. Our model predicts approximately 2- 
3% diffractive events without phase-space cuts, e.g. in diffractive Z-production in pp at ^/s = 1.8 
TeV we obtain 2.64%, where data from DO implies approximately 1.44% m- Restricting the 
phase space in the event generation, by applying the cuts used in the experiments, further reduces 
the fraction of diffractive events, bringing our model closer to data. But the fraction of diffractive 
events is not the complete story. The model should also be able to describe particle spectra, and 
it is thus important to compare the kinematical distributions to data. 


Mean £7- of jets with surviving antiproton Mean tj of jets with surviving antiproton 




Figure 1: Kinematical distributions of diffractive dijet (SD) events compared to non-diffractive 
dijet (ND) events in pp collisions at ^/s = 1.8 TeV obtained with Pythia 8. 

In Fig. [^we show some preliminary results obtained with the new model. We study diffrac¬ 
tive dijet production at the Tevatron, pp —)• pX , [X —)■ JJX '] at y/s = 1.8 TeV. We show the 
mean Et and rj distributions, where the data obtained at the Tevatron showed significant differ¬ 
ences compared to non-diffractive dijet events. SD data revealed a faster falloff in the mean Et 
distribution compared to ND events, and the events were shifted towards positive rj, the proton 
direction. These differences implied a steeper x dependence in the SD events than in ND events. 
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Unfortunately, our model does not capture all of these effects. The SD events generated with the 
new model are boosted towards positive r] which is hne, but the falloff in is not significantly 
steeper than the ND distribution. Our model simply allows for too many high-pT events. While it 
may not solve all problems, we intend to develop a new description of the Pomeron flux to improve 
this spectrum. This should also improve the soft diffraction model implemented in Pythia 8. 

4 Conclusion 

We have presented a review of the soft diffraction models implemented in the general-purpose 
event generator Pythia 8. This soft diffraction machinery allows for QCD interactions and gives 
an decent description of diffractive phenomena. Comparisons to data shows that there is room 
for improvement in the default settings and a new parametrisations of the Pomeron flux is called 
for. A new model for hard diffraction has also been presented, now for the first time allowing for 
non-QCD processes as well has very high-pT’ QCD processes in diffractive systems. The model 
is successful in describing the RGSP, and diffractive fractions obtained with the model agrees 
reasonably with data. Particle spectra obtained with the model has been compared to the data, 
unfortunately not capturing all aspects of the data. Hence the required improvements and updates 
needed in the soft diffraction regime is also needed in the model for hard diffraction. 
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